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Abstract. Time delay measurements have recently been 
reported for the lensed quasar PG 1115+080. These mea- 
surements can be used to derive Hq , but only if we can con- 
strain the lensing potential. We have applied a recently de- 
veloped deconvolution technique to analyse sub-arcsecond 
/ band images of PG 1115+080, obtained at the Nordic 
Optical Telescope (NOT) and the Canada France Hawaii 
Telescope (CFHT). The high performance of the deconvo- 
lution code allows us to derive precise positions and mag- 
nitudes for the four lensed images of the quasar, as well as 
for the lensing galaxy. The new measurement of the galaxy 
position improves its precision by a factor of 3 and thus 
strengthens the constraints on the lensing potential. With 
the new data, a range of models incorporating some of 
the plausible systematic uncertainties yields H — 53ly° 
kms" 1 Mpc -1 . 
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1. Introduction 

Gravitationally lensed quasars with multiple images pro- 
vide a unique tool for determining the cosmological pa- 
rameters independently of the classical methods. In par- 
ticular, time delays between the images of a lensed quasar 
can be used to determine the Hubble constant H (e.g. 
Refsdal 1964, Grogin & Narayan 1996). 
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Extensive photometric monitoring of lensed quasars 
has been undertaken to obtain accurate light curves for 
the most promising objects. Unfortunately, precise time 
delays have been difficult to obtain due to poor temporal 
sampling of the light curves, possible contamination by 
microlensing and short time scale events. 

Light curves have recently been obtained for the 
quadruply imaged quasar PG 1115+080 (Schechter et al. 
1997). The good time coverage, together with the fact 
that microlensing does not seem to significantly contam- 
inate the intrinsic photometric variations of the lensed 
source, yield measurements of two time delays between 
three components of the system (Schechter et al. 1997, 
Bar-Kana 1997). Keeton and Kochanek (1997, hereafter 
KK97) showed that it was difficult to constrain Ho be- 
cause of degeneracies in the lens models. However, KK97 
also showed that the constraints could be improved by 
measuring the position of the lensing galaxy relative to 
the quasar images better than was possible with the pre- 
rcfurbishmcnt HST observations of Kristian et al. (1993, 
hereafter K93). The aim of the present paper is to provide 
a better determination of the galaxy position and, thus, 
an improved estimate of Hq. 



2. Observations 

The observations consist of five I band images of PG 
1115+080 with a typical seeing of / /6-0 / /7. Using SIS, the 
tip-tilt adaptive optics camera installed at the Cassegrain 
focus of the CFHT, two 600s exposures were taken on 
the night of 1995 December 27. The detector was a Loral 

3, 2048x2048 pixel CCD (pixel scale 0'.'0865). In addi- 
tion, three / band frames were obtained at the NOT in 
direct imaging mode. Two 200s exposures were taken on 
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Fig. 1. Top panels: I band images of PG 1115+080 obtained 
at the CFHT (600s exposure, field of 5'.'5) and at the NOT 
(200s exposure, field of 6'.'8). Middle panels: the same images, 
but deconvolved using our new algorithm. Bottom panels: the 
deconvolved lensing galaxy alone. Note that the greyscale was 
chosen to display the full dynamic range of the images. 



the night of 1996 June 7 and an equivalent exposure was 
obtained on 1996 June 10. The Brocam2 2048x2048 CCD 
camera was used (pixel scale 0'.'1071). 



3. Deconvolution - a new algorithm 

The images were deconvolved using a new method, 
described by Magain, Courbin & Sohy (1997), that al- 
lows precise photometric and astrometric measurements of 
strongly blended point sources superimposed on a diffuse 
background. The main idea of this algorithm is to decon- 
volve by a PSF narrower than the total observed PSF so as 
to preserve a good sampling of the deconvolved image. As 
a consequence, 1) the image can be decomposed in a sum 
of point sources, plus diffuse background and 2) the final 
resolution of the deconvolved images is chosen by the user. 



We chose the final deconvolved PSF of the CFHT as 
a Gaussian with a Full- Width-Half-Maximum (FWHM) 
of 3 pixels. The NOT images were deconvolved on a grid 
of pixels two times smaller than the original frames (i.e. 
0'.'05355) in order to keep a good sampling of the lensing 
galaxy, and the final PSF was a Gaussian with a FWHM 
of 6 (smaller) pixels. This led to resolutions of 0'.'26 for 
the CFHT and 0"32 for the NOT. The PSF used to reach 
these resolutions was constructed from two stars about one 
magnitude brighter than the quasar image and situated 1' 
away from the lens. 

The weight attributed to the local smoothing of the 
background component in our deconvolved model image 
(see Magain et al. (1997) for more details) was chosen so 
that the residuals at the location of the lensing galaxy 
had the correct statistical distribution, i.e. Gaussian with 
a standard deviation of 1 (in units of the photon noise) . In 
this way we avoided under- or over-fitting in the area of in- 
terest. A constant smoothing term was used for the whole 
image and was adapted to the area of the lensing galaxy. 
This leads to slight overfitting of the sky noise farther away 
from the target and produces some noise amplification at 
low light levels. 

In the CFHT images, the PSF shows variations of the 
order of 5-10% from one side of the total field to the other. 
Therefore, we allowed the PSF to depart slightly from that 
derived from the two neighboring stars. The PSF was con- 
strained to deviate by the smallest amount consistent with 
an artefact-free deconvolution of the 4 quasar images. 

Fig. 1 shows the result of the deconvolution for a 
CFHT and a NOT image. Since the greyscale was chosen 
to display the full dynamic range, this image clearly shows 
that there are no "deconvolution artefacts" or "ringing ef- 
fects" around the point sources. 

4. Geometry of PG 1115+080 

Five images were used to measure the position of the lens 
galaxy. For each frame, the deconvolution procedure re- 
turned the coordinates of the four point sources, the de- 
convolved image as well as the image of the diffuse objects 
(Fig. 1). The latter was used to measure the position of 
the lensing galaxy. 

To compare our results with K93's results, we trans- 
formed our coordinates to the coordinate system they 
used. The scaling factor and rotation angle were chosen 
to match the positions of the four point sources as well 
as possible, and the transformation was then applied to 
the galaxy coordinates. Fig. 2 compares the new posi- 
tions with the positions from K93. The precision of the 
galaxy position has been increased by more than a factor 
of 3, from 50 mas with the pre-refurbishment HST images 
(K93), to 15 mas with the new data and deconvolution 
technique. Note that K93 used a pixel scale of 0'.'04389, 
which Gould & Yanny (1994) revised to 0'.'04374. Table 
1 summarizes the geometry of PG 1115+080, in the same 
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Fig. 2. Positions of components Al, A2, B and the lensing 
galaxy, relative to component C. The scale for the point sources 
is 10 times larger than for the galaxy. Triangles and squares 
indicate individual measurements from the NOT and CFHT 
images, respectively. Filled circles with error bars show our fi- 
nal results from the weighted mean of the individual measure- 
ments. Open circles with error bars show the HST positions 
from K93. 



orientation as K93, but using the pixel size given by Gould 
& Yanny (1994). 

The galaxy position was measured on each deconvolved 
background by the first two authors independently, us- 
ing either Gaussian fitting or the first order moments of 
the light distribution, resulting in 10 measurements (one 
for each author and each frame). We first averaged the 
two measurements for each frame and then computed a 
weighted mean of the 5 independent measurements, with 
a weight twice as large for the NOT images as for the 
CFHT images because the latter are corrupted by PSF 
variation across the field. 

The standard deviation of the mean is ~ 2 mas for the 
point sources and ~ 10 mas for the lensing galaxy. The ef- 
fect of a variable PSF was tested by performing numerical 
simulations which indicate that possible systematic effects 
could affect the point sources by another 2 mas and the 
lensing galaxy by an additional 10 mas. This leads to final 
lcr error bars of 3 mas for the point sources and 15 mas 
for the lensing galaxy. 

The intensity ratio of A2 relative to Al is given in 
Table 1. This ratio does not show significant variations 
among our frames or between our frames and the K93 
frames, even though they were taken at four different 
epochs. This suggests that microlensing does not signif- 
icantly affect the light curves of Al and A2 in the I band. 



Table 1. Relative positions of the four lensed images of 
PG 1115+080 relative to component C (in arcseconds). The 
intensity ratios in the I band are given for the NOT observa- 
tions obtained on the night of 1996 June 7. The magnitude of 
Al is I = 16.34 ± 0.05 



Image 


I/Iai 


xo 


Y (") 


Al 


1.000 + 0.000 


+1.291 + 0.003 


+2.028 + 0.003 


A2 


0.683 + 0.010 


+1.445 ± 0.003 


+1.578 + 0.003 


B 


0.164 + 0.003 


-0.364 + 0.003 


+1.940 + 0.003 


C 


0.255 + 0.005 


+0.000 ± 0.000 


+0.000 ± 0.000 


Lens 


I = 19.6 + 0.1 


+0.332 + 0.015 


+1.339 + 0.015 



The galaxy magnitude was derived on the NOT de- 
convolved frames by aperture photometry (079 diaphragm 
in diameter). The zero point was computed using several 
standard stars. Images with a much higher signal-to-noise 
would be necessary to derive the shape of the galaxy. How- 
ever, the galaxy observed in the present data is compatible 
with a fuzzy circular object, broader than a point source. 



5. Implications for H 

Schechter et al. (1997) recently obtained light curves for 
images B, C, and the mean of the close pair A=A1+A2 
(Fig. 1). Bar-Kana (1997) analysed the light curves using 
a variety of assumptions about how to treat the photo- 
metric errors. He found that the best resolved time delay 
is Atbc = 25.01"J'y days, and that the other independent 
time delay is best expressed in terms of the time delay 
ratio tabc = Atac/Atba = 1.13i" "_ 17 , with a ~ 0.2 sys- 
tematic uncertainty. 

Using the quasar image and galaxy data of K93, Kee- 
ton, Kochanek & Seljak (1997), Schechter et al. (1997), 
and KK97 found that the system can be fit only by includ- 
ing a perturbation from a small group of galaxies near the 
lens galaxy (Young et al. 1981, Kundic et al. 1997). There 
are significant degeneracies in the models related to the 
position of the group and to the profiles of the galaxy and 
group, but with more precise data the degeneracies can be 
reduced (KK97). In particular, reducing the uncertainties 
in the lensing galaxy position can constrain the position 
of the group. 

We studied the effects of the improved measurement 
of the galaxy position by recomputing the models from 
KK97. We modeled the galaxy as an ellipsoidal mass dis- 
tribution with a variable position, ellipticity, and orienta- 
tion, and we used both an isothermal (dark matter) model 
and a constant M/L model. The constant M/L model 
was approximated by a modified Hubble model because 
it has an analytic deflection formula and is therefore sim- 
pler than a de Vaucouleurs r 1 ' 4 model; the Hubble model 
was chosen to have a fixed core radius s = 0'.'2 (0.55ft. -1 
kpc). Wc modeled the group as a spherical mass distri- 
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Fig. 3. Normalized probability distributions for Ho computed 
from a Bayesian analysis with the four classes of lens models, 
and a total probability distribution taking into account the rel- 
ative probabilities. In the top panel the Bayesian analysis does 
not include the constraint from the time delay ratio vabc, and 
in the bottom panel it does use vabc- The relative probabili- 
ties of the four classes of models, and their implied values for 
Ho, are given in the key. This Figure is to be compared with 
Fig. 6 of KK97. 



bution with a variable position, and we considered both 
a singular isothermal sphere (SIS) model expected for a 
dark matter halo, and a point mass model to examine the 
effects of making the group more concentrated. Increas- 
ing the core radius of the galaxy or group would decrease 
the implied value for Hq, but dark matter distributions 
generally appear to be singular. 

Table 2 compares the new results with the results from 
KK97 to show the effects of improving the galaxy position. 
For models with an isothermal galaxy the best-fit \ 2 has 
decreased despite the fact that the constraints have be- 
come stronger, while for models with a Hubble galaxy the 
X 2 has increased slightly but the fit is still good. 

Following KK97, we determined Hq from a Bayesian 
analysis of the lens models. Using discrete model classes 
limits the generality of the analysis, but until even bet- 
ter constraints are available (such as the shape and profile 
of the galaxy) a full model survey is unwarranted. The 
four model classes examined here illustrate the range of 
effects from model uncertainties; for other model classes 
see KK97. By reducing the range of acceptable lens mod- 
els, the improved constraint from the galaxy position de- 
creases the H error bars by a factor of two (see Table 2 
and Figure 3). The new position also increases the prob- 
ability of the isothermal (dark matter) galaxy relative to 
the constant M/L galaxy. Finally, it slightly strengthens 



Table 2. Model results, including the absolute x 2 °f the 
best-fit model (with Ndof = 3) and the values for Ho from 
the Bayesian analysis including the constraint from the time 
delay ratio vabc- Quantities in parentheses are results from 
KK97 using the data from K93. 

Galaxy: Isothermal Modified Hubble 



SIS Group 


x 2 


1.35 (1.77) 


2.09 (1.47) 




Ho 


49± 7 6 (44 ±11) 


63± 7 6 (55 ± 14) 


Point Mass 


r 


1.32 (1.83) 


1.95 (1.70) 


Group 


Ho 


54±£ (51 ±11) 


73+* (64+1?) 



the discrepancy noted by KK97 between the lens models 
and Bar-Kana's (1997) value for tabc- Improved measure- 
ments of tabc would reveal whether the apparent discrep- 
ancy is due to systematic effects in the measurements or 
to a poor lens model. 

Based on K93 galaxy position, KK97 found Hq = 



511^ kms _1 Mpc _1 . The Bayesian analysis using our 
new data gives a total probability distribution with H = 
53ly kms -1 Mpc -1 , where these error bars incorporate 
the time delay uncertainties as well as some systematic 
uncertainties in the lens models related to the profiles of 
the galaxy and group and to the position of the group. 

Our new deconvolution technique could be used to 
analyse all the images used for the light curves and to 
derive the independent light curves of Al and A2. The 
high precision of the method might, in addition, make it 
possible to further narrow the uncertainties on the time 
delays and hence on H . 
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